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ABSTRACT 

The double pulsar system PSR J0737— 3039A/B is a double neutron star binary, with a 2.4 hr orbital 
period, which has allowed measurement of relativistic orbital perturbations to high precision. The 
low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, 
point to a different evolutionary scenario compared to most other known double neutron star systems. 
We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting 
constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737— 3039A, 
to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 
90 ± 11 ± 5°. Furthermore, we find a mean 95% upper limit on the misalignment between its spin 
and orbital angular momentum axes of 3.2°, assuming that the observed emission comes from both 
magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the 
second pulsar was relatively symmetric, possibly involving electron-capture onto an 0-Ne-Mg core. 
Subject headings: binaries: general — pulsars: general — pulsars: individual (PSR J0737— 3039A/B) 
— stars: evolution 



1. INTRODUCTION 

Pulsars are the rapidly rotating neutron star (NS) rem- 
nants of supernova (SN) explosions, which emit beamed 
and highly coherent radiation, most probably from field 
lines emanating from the magnetic polar regions. This 
radiation can be viewed from Earth as a radio pulse each 
time the beam sweeps across our line of sight. This pulse 
represents the variation in fiux across the portion of the 
emission beam that we observe. 

In general, NSs that are members of binary systems 
represent a separate population from the "standard" iso- 
lated pulsars, which generally have ^ 1 s rotation peri- 
ods. Most pulsars observed in binary systems are "recy- 
cled" to have faster rotation speeds in a phase of mass 
and angular-momentum transfer from the companion 
star. Many neutron star-white dwarf (NS-WD) systems, 
and especially double neutron star (DNS) systems, are 
in close orbits, often displaying relativistic effects. Mea- 
suring the pulse times-of-arrival from these systems has 
provided some of the most stringent and varied tests of 
general relativity (GR) available, as well as the most pre- 
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cise measurements of compact object masses. 

The double pulsar system, PSR J0737-3039A/B, was 
discovered as the most relativistic DNS binary yet to 
be found (Burgay et al. 2003) and, as its nickname sug- 
gests, the only known s ystem in which b oth NSs have 
been visible as pulsars |Lvne et al.l l2004f ) : one is a re- 
cycled pulsar, PSR J0737-3039A (henceforth referred 
to as pulsar "A"), with a spin period of 22.7 ms; the 
other pulsar, PSR J0737-3039B (henceforth referred to 
as pulsar "B"), which has a 2.77-s rotation period, was 
formed in the second SN and has not been spun up. The 
highly relativistic nature of this system, and the pres- 
ence of two pulsars, has provided pulsar astronomers 
and NS theorists with many astrophysical phenomena 
to study in greater detail, and with more precision, than 
ever before: it immediately provided the b est prediction 
of the coalescence rate of DNS systems (jBurgav et al.l 
120031: iKalogera et al.|[200l : it also allowed for the mea- 
surement of the spin-orbit precession of pulsar B by 
modeling the eclipse light curves as pulsar A passed be- 
hind pulsar B, con sistent with the predictions of GR 
(jBreton et al.l 12008"). A more recent study of pulsar B 
has constrained its pulse profile evolution due to geodetic 
precession, thus constraining models of its emission beam 
shape, and has confirmed that, due to precession, our 
line of sight cu rrently int ersects an em pty portion of the 
beamJPerera et al.ll2010f l. In contrast. [Manchester et al.l 
ll2005i) found no evidence for secular change in the pulse 
profile of pulsar A over a 3 year data span, suggesting ei- 
ther a close alignment of the pulsar spin axis with the 
orbital angular momentum, or a fortuitous precession 
phase. 

Perhaps the best-known result from observations of 
the double pulsar is the long-term timing analysis 
that has led to the confirmation of the predictions of 
GR in the strong-field regime to within 0.05%, based 
on the measuremen t of Shapiro delay in this system 
([Kramer et al.l l2006f l . The parameters describing the 
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PSR J0737— 3039A/B system continue to be monitored 
and updated, providing increasingly rigid constraints on 
relativistic theory. An extensive review and specula- 
tio ns on future science wi th the double pulsar is given 
bv lKramer fc StairsI (l200l . 

1.1. A Different Evolution? 

Several of the properties of the PSR J0737-3039A/B 
system differ significantly from other known DNS bi- 
naries. In particular, the low pulsar-B mass of 
1.2489(7) Mq, the small measured eccentr icity (0.0878) 
and low transverse velocity (lOkms"^) (jKramer et al.l 
l2006f ). show a marked departure from some DNS sys- 
tems for which these parameters have been measured, 
though also similarities with others, as discussed be- 
low. This is certainly due to the specific evolu- 
tion of these systems. It is generally thought that 
DNS systems result from one of three evolutionary 
scenarios, which we briefly describe here. A more 
detailed overview of DNS an d binary NS evolution 
can b e found in, e.g., Bhattac harva fc van den Heuvell 
([T99lh. IPh inncv fc Kulkarni f 19941) . IStairsI (|2004D . and 
iTauris fc v an den Hcuvcl (2006). 

In the "standard" scenario, the system initially con- 
sists of two massive stars, typically > 8 — ^Mq. The 
more massive of these eventually evolves off the main 
sequence, at which time it fills its Roche lobe, stably 
transferring matter to the secondary star, before col- 
lapsing to explode in a SN, leaving behind the first- 
formed NS in a high-mass X-ray binary (HMXB; see, e.g., 
ITauris fc van den Heuvell 120061) . The remaining massive 
secondary then evolves to overfill its Roche lobe. Un- 
stable mass transfer then occurs, creating a common en- 
velope (CE) that engulfs the NS, which then spirals in 
toward the companion, picking up angular momentum 
and spinning up as a result. Dynamical friction results 
in the ejection of the CE, leaving behind a NS with the 
remaining He core of the companion star. Further mass 
transfer may ensue, in which case there is a further in- 
crease in the NS rotation speed. The He star then even- 
tually undergoes a SN, leaving behind the DNS. 

The double-helium core scenario was hypothesized to 
avoid possible hypercritical accretion onto the NS during 
the CE phase that may cause colla pse into a bla ck hole, 
rendering the system unobservable (|Brownlll995[) . Here, 
the two massive progenitor stars initially have masses 
within ^ 10% of each other, and the orbit is sufficiently 
wide for the primary to evolve a CO core. The secondary 
then evolves off the main sequence before the primary 
undergoes a SN. The mass transfer rate in this case is 
very high, such that the CE that develops is comprised 
of both stars' envelopes. This causes a rapid spiral-in of 
the cores and CE ejection. After the primary undergoes 
SN, the remaining evolution resulting in the DNS is sim- 
ilar to the standard process described above. Although 
it is expected to occur 2 — 10 times less f requently than 
in the standard picture (jDewi "eralll2006l) . it may be the 
predominant formation channel; this depends on the like- 
lihood of the standard scenario leading to black hole for- 
mation due to hypercritical accretion by the fi rst-formed 
NS in a CE, thus inhibiting DNS formation ([Chevalieii 
119931) . On the other hand, it has been shown in the 
case of the NS-CO white dwarf binary PSR J1802-2124 
that it is possible for a close NS binary system to sur- 



vive a CE phase without hypercritical accretion occur- 
ring (jFerdman et al.|[2010l and references therein). 

A third scenario, commonly referred to as an electron- 
capture supernova (ECS), separates itself in how the sec- 
ondary star evolves to form a NS. Here, an 0-Ne-Mg 
core passes a threshold density that allows electrons to 
be captured on ^^Mg. This decreases the electron degen- 
eracy pressure, which in turn lowers the Chandrasekhar 
mass of the core, inducing its collapse (iMiyaii et al.llT980l : 
iNomotol 119841 : iPodsiadlowski et al.l [20051). This would 
result in low NS velocities, since this type of event 
is thought to proceed over a much shorter timescalc 
than that needed to develop the instabilities thought 
to cause substantial S N kicks (jPodsiadlowski et"ani2004l : 
Ivan den Heuvell [20041 ) . and has also been invoked to ex- 
plain the formation of a subset of NSs with significantly 
lower measured masses relative to the overall known 
NS mass distribution ( including PSR J073 7-3039B, dis- 
cussed in this work; [Schwab et all [2010[ ). The ECS 
scenario is also a suspected mechanism for forma- 
tion of wide-orbit H MXBs consisting of Be-star donors 
([Linden et al.ll2009D . 

In order to investigate how these scenarios may apply 
to DNS systems we observe, we list those with measured 
masses in Table [TJ and the values of several of their prop- 
erties. ^PSR^^2127-fllC resides in the globular cluster 
Ml 5 ([Jacobv et al.l I2006I ): the evolution of this pulsar 
thus likely involved unique dynamics such as exchange 
interactions that are extremely rare even in the densest 
areas of the field. It is thus reasonable to assume that the 
evolution of this pulsar should differ substantially from 
the others in Table [TJ PSR J1518+4904 is also a curious 
system, in that it has a moderate eccentricity, compara- 
ble to that of PSR B1534-I-12, however, its measured 
transverse velocity is relatively low. Most interesting 
is that the component masses seem to be relatively ex- 
treme, with a very light first-formed NS and a relatively 
heavy companion NS. This, along with the system's rela- 
tively long orbital period, indicates that it may have pro- 
ceeded through an entirely unique e volutionary scenario , 
which is not yet well understood ([Janssen etaDHool 
[Wong et al.|[20T0l) . PSR J1829-h2456 may be a close 
cousin to PSR J1518-f4904, given its currently known 
mass limits, as well as its spin and eccentricity, which 
are s imilar to the latter system ([Champion et al.l [2004 
120051 ). PSR J1811-1736 has the most eccentric orbit in 
Table [U indicating it may have narrowly escaped disrup- 
tion after the second SN event. It is unlikely that reliable 
proper motion estimates will be obtained in the near fu- 
ture, given the large estimated distance to this pulsar, 
along with its relatively poor timing precision. Such a 
measure ment would help to better discern its formation 
history (|Lvne et al.ll200l3HCorongiu et al.ll2007D . Finally, 
PSR J1906-h0746 was originally thought to be a DNS 
system in which we see the low-m ass young companion 
NS to an un seen recycled pulsar ([Lorimer et al. I [20061 : 
lKasian|[2OO80 . The pulsar in this system is believed to 
be a young NS. Following updated mass measurements 
from longer-term timing analysis (Kasian 2012), the na- 
ture of the companion is now ambiguous, and it may well 
be a massive WD. In this case, it may have undergone 
an evolution similar to that of PSR J1141— 5645, another 
close NS-WD binary in a mildly eccentric or bit in which 
the massive WD was formed before the NS ([Kaspi et al.l 
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TABLE 1 

Known Double Neutron Star Binary Systems with Measured Component Masses. 



Pulsar 


p . 

^ spin 


e 


^orb 






Vtr 


P b 


References 




(ms) 




(days) 


(Me) 


(Mq) 


(kms-1) 


(years) 




J0737-3039A/B<: 


22.70/2773 


0.0878 


0.102 


1.3381(7) 


1.2489(7) 


10<^ 


75 


1 


J1518+4904 


40.94 


0.249 


8.63 


< 1.17 


> 1.55 


25 




2 


B1534+12 


37.90 


0.274 


0.421 


1.3332(10) 


1.3452(10) 


107 


706 


3 


J1756-2251 


28.46 


0.181 


0.320 


1.312(17) 


1.258(18) 


C 


488 


4 


J1811-1736 


104.2 


0.828 


18.8 


< 1.64 


> 0.93 






5,6 


J1829+2456 


41.01 


0.139 


1.18 


< 1.34 


> 1.26 






7,8 


J1906+0746f 


144.1 


0.0853 


0.166 


1.290(11) 


1.323(11) 




169 


9,10 


B1913+16 


59.03 


0.617 


0.323 


1.4408(3) 


1.3873(3) 


88 


296 


11 


B2127+11C 


30.53 


0.681 


0.335 


1.358(10) 


1.354(10) 


168 


278 


12 



Ref erences. — (1) IKramer et al.l |(2009); (2) IJanssen et al.| l(2003); (3) IStai rs et a ll ll200 l): (A) IFerdmanI 
1I2OO8I'): (5UL vnc ct al. (2000); (6) Corongiu ct al. (2007); (7) Champion et al. (2004); (8) Champi on"et al.l (120051 ): 
(9) ILorimer e t al. (2006); (10) Kasian (2012); (11) Weisberg et al. (2010); (12) Jacoby et al. (2006). 

Here, mi and m2 refer to the first-formed (usually recycled) and second-formed (usually unrocycled) neutron stars, 
respectively. In all cases, the first-formed NS is observed as a pulsar, except in the case of PSR J0737— 3039A/B, 
where both have been seen as pulsars. 

Precession period of the first-formed neutron star in the system. 

For the double pulsar we show values for both pulsars where they differ. 

This value is based on a distance of 500 pc to the pulsar, based on the dispersion measure found from timing 
analysis. 

So far only proper motion in the RA direction has been measured to be —0.7(2) masyr~^, corresponding to a 
velocity in the RA direction of ~ 8.3kms~^. 

' It is now thought that the companion to PSR J1906-I-0746 is likely to be a high-mass WD, however, a NS companion 
is not ruled out. 



[2OOOI: iTauris fc Sennelsl [2Q00l: IBhat etHI 120081) . It has 
not yet been ruled out, however, that this system is a 
DNS. If this is the case, an ECS may be the prevalent 
scenario for formation of the first-born NS; an ECS is 
also a potential channel for the second SN given the rel- 
atively low system eccentricity, alt hough the co mpanion 
mass may exclude this possibility ()Kasianll2012D . 

One can see similarities between certain systems. In 
particular, PSR B1913+16 and PSR B1534-1-12 have 
massive companions, large eccentricities, and high trans- 
verse velocities. It is believed that these systems have un- 
dergone either a standard or double-He core-collapse SN 
of an iron core in forming the second NS. Such high eccen- 
tricities and space velocities are indicative of a high mass- 
loss, asymmetric SN from a massive p rogenitor, giving 
a significant iiatal k i ck to the system (iWex et al.ll200Cll: 
Willems et all [200I iThorsett et all l2005t iWong et alJ 
20!^^^ 

In contrast, the PSR J0737— 3039A/B system has an 
order-of-magnitude smaller eccentricity and transverse 
velocity than those measured for t he PSR B1 5 34-I-12 
and PSR B19134- 16 systems (jStairs et al.l 120021 : 
iWeisberg et al.l l2010t ). This indicates there may have 
been little mass lost, and a small kick to the system 
from the second SN. In addition, the mass of pulsar B 
is also significantly lower than any of the components 
in these two other DNS binaries. Taken together, these 
properties suggest a different evolutionary path has been 
taken in the formation of pulsar B. 

It has been suggested that this may best be ex- 
plained by pulsar B hay ing formed after an ECS 
(jPodsi adlowsk i et al.l I2005D . In addition to the low 
transverse velocity of the PSR J0737-3039A/B sys- 
tem, the mass of pulsar B provides important evidence 
for this. The mass corresponding to the binding en- 
ergy of the NS for many equati ons of state is given by 
~ QMAjMnf^/Mp^Y M p,c^ (|Lattimer fc Yahill Il989l: 
iLattimer fc PrakashI 120011) . corresponding to a mass 



Mbe ^ 0.13 A/q for pulsar B (see also IStairs et alJ 
|2006[ ). Through independent modeling of the pre-SN 
mass of the B p ulsar based on a collapsing 0-Ne-Mg core, 
Podsiadlo wski et al.l (|2005l ) have found that the critical 
mass for collapse should range from 1.366 to 1.375 M©. 
Compa ring with the measu red mass of 1.249 Mq for pul- 
sar B (jKramer et al.l 120061 ) . this is consistent with the 
above estimate for Mbe- This would mean that almost 
no baryonic matter has been lost during the SN event. 
Any remaining energy would have gone into changing 
the orbital properties of the system, or contributing to a 
kick at the time of the SN. This is not only evidenced by 
the low measured projected velocity of the system, but 
also its small eccentricity. 

Finally, we leave discussion of the PSR J1756--2251 
system, and how it compares to the double pulsar sys- 
tem, until SectionlHl However, one can see that the com- 
ponent masses, as well as the orbital eccentricities, are 
very similar to those in PSR J0737— 3039A/B. It can thus 
be argued that this system may have followed a similar 
evolutionary path to that of the double pulsar system. 
This would have important consequences for population 
studies of DNS systems. 

2. UNIQUE EVIDENCE FOR A SYMMETRIC SUPERNOVA 

Together, the properties of the PSR J0737-3039B sys- 
tem present tantalizing clues that the progenitor of the 
B pulsar may have ended its life in an ECS, or other 
symmetric SN event. However, it may still be possible 
to devise an evolutionary model that can reproduce its 
measured properties, especially given the possibility that 
a small transverse velocity does not necessarily imply a 
small velocity in the radial direction, and thus a small 
overall space velocity (jKalogera et al.ll2"008| ). As we now 
discuss, measurement of the spin direction of pulsar A 
can help to decide on the validity of an ECS or iron core- 
collapse scenario as the pulsar B formation mechanism. 

It is expected from accretion theory that the spin 
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axis of pulsar A will become aligned with the to- 
tal angular momentum of the binary system (well ap- 
proximated by the orbital angular momentum) as it 
accreted matter from the progenitor of pulsar B. If 
the second SN is close to symme tric, this alignment 
would not be disturbed tPodsiadlowski et alJi200 4'). By 
contrast, if there is a large kick to the system, the 
resulting misalignment would equal the angular dif- 
ference in orbital pl ane orientation before and after 
the SN event (e.g., 'Wex et al.' '2000"). Several stud- 
ies (e.g.. [P ewi fc va n den Heuvel 200 4; Piran fc Shav i^ 
200l iWihems et al.ll2006t I'Stairs et'ani2(MIWong et all 



20101) have examined the explosion of the pulsar B pro- 
genitor. Using the t iming-derived prop er motion from 
iKramer et all (|2006[ ). [Stairs et al.l (|2OO60 predict a post- 
SN m isahgnment an g le S f or PSR J0737-3039A of < 
11°. Willcms et al.l ([2005') predict compatible values 
under different kinematic and progenitor-mass assump- 
tions. 

If the above studies are correct (see iKalogera et al.l 
120081 for a discussion and criticis m of the assuniption 
of a small radi al velocity used in iPiran fc Shavivll2005l 
and iStairs et a l. 2006), the measurement of a low spin- 
orbit misalignment angle for PSR J0737— 3039A — in con- 
junction with the low pulsar B mass, system eccentric- 
ity, and transverse velocity that have been determined 
from timing measurements — would thus provide very fa- 
vorable evidence for an ECS/symmetric SN event having 
formed pulsar B. We thus aim to determine the orbital 
geometry of the PSR J0737-3039A/B binary in order to 
constrain the misalignment between its orbital angular 
momentum and the pulsar A rotation axes. 

One way to accomplish this is through investigation of 
the observable effects of geodetic precession on pulsar A. 
In general, precession of the spin axis orientation of a 
pulsar causes a change in our line of sight through the 
emission region over time. This is reflected in a contin- 
uous change in the observed pulse profile shape. The 
extent to which these effects can be observed depends 
on the spin and orbital geometries of the pulsar system. 
This is particularly the case for the misalignment angle 
S between the pulsar spin and orbital angular momen- 
tum, which forms the opening angle of the cone that is 
swept out by the spin axis over the course of a preces- 
sion period. These geometric parameters can be modeled 
and determined through long-term analysis of the pulse 
profile. 

The precession periods for the other relativistic 
DNS binaries mentioned earlier, PSR B1534-I-12 and 
PSR B1913-f 16, are approximately 706 and 298 years, 
respectively. Misalignment angles have been measured 
for these pulsars: S = 25.0 ± 3.8° and 18 ± 6° 
for PSR B1534+12 and PSR B19 1 3-H6, respectivel y 
(jWeisberg et al.l [19891: IkTm^ [19981 : TStairs et al.ll27)0l . 
By comparison, PSR J0737— 3039 A has a precession pe- 
riod of 75 years (|Lvne et al.ll200l . The effects of geode- 
tic precession on the pulse profile in this system should 
thus be more readily observable on a shorter timescale 
than for PSR B1913-I-16 and PSR B1534+12, provided 
that S is non-negligible. 

In this article, we present a study of profile evolution 
for pulsar A over 6 years of observation. After giving a 
brief overview of previous work in determining the geom- 
etry of pulsar A in Section [3l we describe our observa- 



tions in Section 131 We outline our analysis and calcula- 
tions of pulsar A profiles in order to obtain geometrical 
constraints in Section [5] In Section |6| we describe our 
results. In Section [8| we discuss the implications of these 
results and compare them to previous findings. Finally, 
in Section [9| we provide concluding remarks. 

3. PREVIOUS WORK 

Radio polariz a tion data analysis was performed by 
iDemorest et al.l (|2004[ ) in order to derive geometrical 
constraints from position angle variations across the 
PSR J0737-3039A profile. That work found solutions 
that prefer a spin-magnetic axis separation angle value 
of a 5 ° or a ^ 90 °. They favor the former low-a so- 
lution, in order to a void an emission beam opening angle 
^ 180°. However. IDemorest et al.l (|2004[ ) describe the 
changing position angle across the profile using the rotat- 
ing vector model (RVM; .Radh akrishnan fc Cookel [l969). 
which does not appear to be a good fit to the data (by 
the authors' own admission, the RVM provides an unsat- 
isfactory descri ption of the profiles o f many recycled pul- 
sars; see, e.g., Navarro et al. 1997|; iXilouris et al.l [19981 : 
[stairs et al., 1999: Yan ct al. 201^. 

Using a data set th at spans almost three years, 
[Manchester et al.l ([2005D found no evidence that the pro- 
file width of PSR J0737-3039A was changing with time. 
Because of the symmetric properties of the observed 
pulse profile, they interpreted both pulse components as 
emission from a single magnetic pole. To avoid implying 
a beam opening angle of greater than 90 °, they suggested 
a small but non-zero spin-orbit misalignment angle, with 
5 ~ 14 ° as the preferred value, combined with a pre- 
cessional phase close to either 0° or 180°. A non-zero 
misalignment angle for pulsar A implies that pulsar B 
received a significant natal kick. However, the subse- 
quent precision timing analysis by IKramer et al.l (|2006l ) 
and very long baseline interfero metry measu r emen t of its 
annual geometric parallax by iDeller et al.l (|2009[ ) have 
both confirmed the proper motion, and thus transverse 
velocity, of the double pulsar system to be relatively low 
compared to oth er DNS systems . Thi s and the stud- 
ies co nduct ed bv IPiran fc Shayivl (p()05h . [wlliems et al.l 
(|2006[) . and IStairs et al.l (|2006D suggest that a large na- 
tal kick is unlikely for the PSR J0737-3039A/B system, 
assuming the system has a small radial velocity. This 
would in turn predict a low value for 6. 

If PSR J0737— 3039A is a near-orthogonal rotator, then 
the profile we observe from this pulsar may represent 
contributions from each of the two magnetic poles. Ad- 
ditionally, the two pulse components are separated by 
approximately 180 ° and the region between is largely 
free of pulsed emission. With this model there is not 
a problem in postulating a small or even zero value of 
the spin-orbit misalignment angle. A detailed study of 
the pulse shape and polarization, and investigation of 
possible emission mechanisms is necessary to achieve a 
more satisfactory understanding of the geometry of the 
system. 

To further in vestigate, we hav e cond ucted a sim- 
ilar analysis to [Manchester et all (|2005l ). performing 
a fit of the geometrical parameters describing the 
PSR J0737-3039A/B system to the pulse width data. 
We have extended our data set to include 6 years of ob- 
servations of this system, making the constraints on the 
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Fig. 1. — Integrated pulse profiles and difference profile residuals over nearly 7 years of GBT observations. Left: high signal-to-noise, 
aligned and averaged profiles are shown from each of 13 concentrated observing campaigns, and labeled with the central date of the given 
epoch. Right: the most recent of these, taken in 2011 June, is subtracted from each of the other averaged profiles. The difference profiles 
in the right-hand panel are plotted at twice the vertical scale of the left-hand panel. 



system geometry far more rigid 0. As we describe be- 
low, we perform this investigation assuming both one- 
and two-pole beam models; as we discuss in Section [51 
we find that the latter emission structure allows for a 
smaller beam size, and is our preferred description for 
this system. 

4. OBSERVATIONS 

Observations of PSR J0737-3039A were taken be- 
tween MJD 53521 (2005 May 31) to MJD 55731 (2011 
June 19) using the 100-m Robert C. Byrd Green Bank 
Telescope (GBT) in West Virginia, with the Green 
Bank Astronom ical Signal Processor backend (GASP; 
IDemorestl[200l . GASP is a flexible baseband system, 
which performs 8-bit Nyquist sampling of the incoming 
data stream at 0.25/xs intervals in both orthogonal po- 
larizations. The i ncoming data stream was then coher- 
ently dedispersed (Han kins fc Ricket^l 19751) in software. 
After this, the signals were folded at the topocentric 
pulse period of PSR J0737— 3039A to form pulse profiles, 
each of which were accumulated over approximately 30 s. 
This choice of integration time reflects a compromise be- 
tween ensuring a fully-sampled orbit for the very precise 
Shapiro delay measurement for timing studies, and ob- 
taining pulse profiles with adequate signal-to-noise ratio 
to effectively carry out that analysis. 



These folded pulse profiles were typically flux- 
calibrated in each polarization using a reference signal 
from a noise diode source that was injected at the re- 
ceiver. When calibration data were not available, we nor- 
malized the profile data in each polarization by the root 
mean square (rms) value of the corresponding off-pulse 
signal. The data were finally summed over both polariza- 
tions and across all frequency channels to give the total 
power signal for each integration (for f urther details on 
GASP operatio n and data reduction, see lDemorestll200"7l : 
lFerdmanl[2?)08h . 

As stated above, the observing strategy for the dou- 
ble pulsar has been such that it would benefit the timing 
analysi s that has produced the most recent stringent test 
of GR (H ramer etall [20061 ). To that end, we have been 
taking data in monthly observing sessions, each consist- 
ing of a 5-8 hr track of the object. The incoming sig- 
nal was divided into a maximum of 16 and 24 x 4 MHz 
channelj^. centered at 820 and 1400 MHz, respectively, 
in alternate months. This was done in order to effec- 
tively constrain astrometric parameters such as position 
and proper motion, as well as possible changes in the 
measured dispersion measure (DM) of the pulsar. We 
have also conducted semi-annual concentrated observing 
campaigns that cover approximately 1-2 weeks each, and 
which typically contain up to six observing sessions that 



^ A preliminar y version o f this analysis can be found in lFerdmanl 
ifSOOl ) and Fc rdman et al.l pOOS ) 



" The number of channels used occasionally varied due to radio 
frequency interference (RFI) and available computing resources. 
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each last 5-8 hr. The purpose of these was to better con- 
strain orbital parameters and to investigate the effects of 
special relativistic aberration on the pulse shape over an 
orbital period. These campaign observations were taken 
exclusively using the 820-MHz receiver. 

The data set used for this work consists solely of these 
concentrated campaign observations. We have excluded 
the P arkes telescope data set used in [Manchester et all 
()2005( ) , as they were taken with a different observing sys- 
tem, at a different frequency. However, given that there 
is also a lack of change in those pulse profiles, the upper 
limits derived for S in this work can be considered to be 
conservative. 

5. LONG-TERM PROFILE ANALYSIS 

In order to perform a long-term profile shape analysis, 
we first combined the integrated, calibrated pulse profiles 
from each set of campaign observations described above 
into high signal-to-noise profiles. Each of these averaged 
profiles represented the midpoint in time during which 
each corresponding set of observations were taken. The 
left panel of Figure [T] shows the combined average pro- 
files from each representative epoch. For illustration, the 
right panel of Figure [T] shows the residual profiles from 
each corresponding epoch, obtained when the profile rep- 
resenting the data taken around MJD 55721 (2011 June 
9) is subtracted from each of the other mean profiles. 
A qualitative glance at these difference profiles shows 
the lack of substantial changes between epochs. Further- 
more, there seems to be no secular change over time, 
which one may otherwise expect for this highly relativis- 
tic system, and which is the case for both PSR B1913-I-16 
and PSR B1534+12. This would support the hypothesis 
of pulsar A's angular momentum being aligned with that 
of the system. 

5.1. Width Calculations 

After fitting for the pulse height with a simple Gaus- 
sian to several points around the profile peak, we calcu- 
lated the widths and width uncertainties of the pulsar A 
profile at each epoch, at several fractional pulse heights. 
We did this through a bootstrap-type method, as fol- 
lows: we first performed a polynomial interpolation to 
find the pulse phases on the profile, at the chosen frac- 
tional pulse height. This was done for each side of the 
profile (or profile component, in the two-pole case) by 
randomly choosing 11 out of 25 possible (and not nec- 
essarily contiguous) data points along the profile, about 
the location in phase it is expected to cross the given 
height. The difference between the fit phase on each side 
of the profile was taken as the calculated pulse width. 
After 32768 such iterations, we constructed a histogram 
of these trial widths, to which we fit a Gaussian profile. 
The mean value and la width of this distribution were 
taken as the final pulse width and uncertainty for the 
given epoch, respectively. 

To include the possibility that we may be seeing two- 
pole emission from this pulsar, we also treated each com- 
ponent, separated from each other by a region of negligi- 
ble pulsar emission, as the distinct emission from each 
magnetic pole, and separately measured the width of 
each of these as described above. 

To ensure that we obtained a reliable fit, we measured 
the profile widths at 5% fractional height increments. 
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Fig. 2.— Width measurements for PSR J0737-3039A at 40% of 
the pulse height. The top panel shows the widths when assuming 
a one-pole emission model for this pulsar. The bottom two panels 
show the measurements in the case of a two-pole model, where each 
profile component is assumed to originate separately from each 
beam. For this case, the middle and bottom panels correspond to 
the first and second components, respectively. 

between 30% and 50% of the total pulse height. Within 
this region, there are no features (e.g. the plateau re- 
gion in the first pulse component at ~ 10% pulse height) 
that may affect our analysis from obtaining a clear de- 
termination of the pulse width. As an example of our 
width measurements. Figure [5] plots those made at 40% 
of the pulse height for each of the one- and two-pole 
emission cases. At a glance, we see no perceivable sec- 
ular change in the pulse shape. Any perceived outlying 
measurements are likely due to several possible effects; 
principal among these are likely low-level radio frequency 
interference (RFI), and lack of accounting for the non- 
orthogonal polarization feeds in our calibration (see, e.g., 
Ivan StratenI l2006l |2012|) . We performed the geometri- 
cal model fit using the set of measurements at each of 
the fractional pulse heights for which we have calculated 
widths. 

5.2. Geometry Fit 

In order to model the geometry of 
PS R J0737-3039A/B, we use the formulation given 
by iRafikov fc Lail ()2006D . which relates the long-term 
profile evolution of pulsar A to the system geometry. 
The dependence of the observed pulse width on the 
system geometry can be expressed by the following: 



cos <i>o 



cos p — cos C cos a 
sin C sin a ' 



(1) 



where $o is the half-pulse width, a is the angle between 
the spin and magnetic axes, is the angle between the 
spin axis and the observer line of sight, and p is the 
half-opening angle of the part of the emission cone cor- 
responding to the given fractional pulse width. In using 
this equation, we are assuming a circula r emission beam. 
In their work, [Manchester et all (|2005[ ) investigated the 
effect of using a noncircular beam on their analysis, and 
found that the results were only marginally affected. We 
thus believe that a circular emission beam is a reasonable 
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assumption for this analysis. For schematic diagrams 
that show the angles involved in pulsar spin and orbital 
geometries in the a bove expressions, refer to Figure 1 in 
iDamour fc favloil ([l99l . 

In relativistic systems, we would expect that geode- 
tic precession would have the general effect of varying 
the angle C over time. In order to express Equation [T] 
in terms of parameters that do not vary with time, we 
make use of the fact that the direction of the spin axis 
vector si can not only be described by the polar angles 
(Ci v) (where 77 is the angle between ascending node and 
the projection of the spin axis on the plane of the sky), 
but alternatively by {S,(pso), assuming that the preces- 
sion phase varies linearly in time. Here, 6 is the spin- 
orbit misalignment angle, and (j^soit) is the longitude, or 
phase, of si in its precession around the orbital angular 

momentum vector fc, a s measured from the — J axis, as 
illustrated in Figure [3] (jPamour fc Tavlori[T99l . 

The transforma t ion ( C, ij) — >■ {S, (j)so) is given by 
IDamour fc Tavloil (|T99l : 



cos A = cos S cos i — sin S sin i cos (/)so 
sin (5 sin (/)so 



COST] 



lA 



where X — n — (^, and so sin A = sin C and cos A = 
We thus have: 

cos C = — cos S cos i + sin 5 sin i cos (f>so 

sin S sin 0so 

cos 77 = — . 

sm^ 

We can also recast (pso as: 

0so = f^f"(^-ri), 



(2) 
(3) 
cos(. 

(4) 
(5) 

(6) 



where f2^^'" is the angular precession frequency, t is the 
epoch, and Ti is the reference crossing time of the spin 
axis through precession phase 0so =0°. This analysis 
was performed separately for the c ases of cost > (i = 
88.7°) and coai < {i = 91.3°) ()Kramer et all [20061 ). 
We use a single value for the inclination in our analysis 
for each of these two cases; the fractional uncertainty 
in these values is small enough that doing so would not 
significantly affect our results. 
Finally, (C, rj) are now given by: 

cos C = — cos S cos i + sin 6 sin i cos [17^'''" — Ti)] (7) 



cos r/ 



sin^ 



(8) 



Substituting Equation [7] into Equation [T] gives us an 
expression for the time- varying pulse width $0, and that 
depends on the three angles a,d, and p, as well as the 
epoch of zero precession phase Ti. 

These parameters can be determined through a least- 
squares fitting of our pulse width data to the model 
given by Equation [TJ Here, we have performed a three- 
dimensional grid-search fit, over the ranges 0° < a < 
180°, < 5 < 90°, and with Ti running over one pre- 
cession period centered on 1990 January 1 (MJD 47892). 
Although it is physically and mathematically possible for 
S to have values between 90 ° and 18 °, we focus on val- 
ues < 90°; it has been shown bv iBailesI (|1988[ ) that. 
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Fig. 3. — Pulsar precession geometry for a general system. The 
vectors / and J define the plane of the sky, so the K points in the 
direction away from the observer, i and j define the plane of the 
pulsar orbit, so that k points in the direction of the orbital angular 
momentum, i is the orbital inclination relative to the observer 
line of sight. Due to geodetic precession, the pulsar spin axis si 
will trace out a cone with opening angle 5, the misalignment angle 
between the pulsar spin axis and the orbital angular m omentum. 
This diagram is adapted from IDamour & Tavlon l| 19921 1. 



based on modeling of PSR B1913-I-16, misalignment an- 
gles with these values are much more physically likely, 
unless the SN kick was extremely large. 

For each fractional pulse height, we run a Levenberg- 
Marquardt fit of the corresponding width data to the 
above model at each grid point, allowing only p to vary. 
In the one-pole model case, there is only one beam and 
thus a single p value found at each grid point. In the 
two-pole case, for each fractional pulse height, we have 
two width data sets that we fit simultaneously into the 
single model described above. We make the assumption 
that the magnetic poles are separated by 180°, so that 
we substituted a 180 ° — a in the model when describ- 
ing the second component in the pulsar A profile. Here, 
there are two emission beams and thus two distinct p val- 
ues for which to fit (henceforth referred to as pi and p2, 
corresponding to the first and second profile components, 
respectively) . 

After performing a fit over a coarsely sampled grid in 
each dimension, we did a final fit with variable grid-point 
bin sizes, such that regions of higher probability density 
were re-sampled with finer bins, and regions with rela- 
tively little or negligible probability density were probed 
with coarser sampling. This was done to provide an accu- 
rate calculation of probability density, while minimizing 
computation time. We then arrived at a joint probabil- 
ity density distribution for a, S, and Ti, from which we 
found individual probability density functions (PDFs) for 
each of these parameters in turn by marginalizing over 
the other two parameters. These fits and PDF calcula- 
tions were done for both the one- and two-pole emission 
models, for each set of fractional pulse height width data, 
and for the two possible inclination values: i < 90 ° and 
i > 90°. In Figure m we show the resulting PDFs for 
the one- and two-pole model fits to the width data taken 
at 40% of the peak pulse height. We found PDFs for p 
by calculating a histogram of all fit p values, weighted 
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Fig. 4. — Probability density functions (PDFs) for each of the geometric parameters in the model of IRafikov fc Lail II2006I ) for the 
PSR J0737— 3039A/B system, found using widths calculated at a 40% fractional pulse height. Results for the case of orbital inclination 
j < 90 ° are represented by solid lines, and those for i > 90° by dotted lines. Shown are the PDFs for the angle cx between the spin and 
magnetic axes of pulsar A (top left), the misalignment angle 5 between the pulsar A spin axis and the total angular momentum of the 
system (bottom left), the epoch T\ of zero geodetic precession phase (top right), and the half-opening angle p of the region of the pulsar A 
beam that corresponds to the measured width at 40% of the total pulse height (bottom right). The results for the one- and two-pole cases 
are shown in the top and bottom sub-plots for each parameter, respectively. The vertical dashed lines in the T\ PDF plots represent the 
range of our data span. 



by the output probability density at the corresponding 
grid point. Results for both orbital inclinations are very 
similar, which is unsurprising, given how close the true 
inclination is to 90°. Even with independent measure- 
ment of any of the geometrical parameters, it would in 
practice be very difficult to discern which inclination is 
correct; such a measurement would have to be extremely 
precise. In what follows, we will refer solely to the results 
obtained for the case of i < 90 °. 

6. RESULTS 

We present the fit geometric parameters in Table [21 
where our results are shown for each of the one- and 
two-pole models, and for each fractional pulse height at 
which we calculated the input width data. As stated 
above, we only show the results for z < 90 °. The median 
and 68% confidence intervals for the histogram are used 



as the measured a and p values and their corresponding 
uncertainties, respectively. For 5, we have determined 
68%, 95%, and 99% upper limits. 

The one- and two-pole emission cases give consistent 
results for a and 5, favoring the pulsar A axis of rotation 
to be at approximately right angles to the magnetic axis, 
and a low spin-orbit misalignment angle. However, a 
major difference between the two emission models comes 
with the preferred values we find for p. Since p {pi and 
P2 in the two-pole case) is the beam half-opening an- 
gle corresponding to the given fractional pulse width, it 
is always smaller than the half-opening angle of the full 
beam. The values we find for p in the one-pole case thus 
greatly favor a full beam opening angle that is larger than 
180°. The two-pole case, however, easily allows for two 
emission beams on opposite sides of the NS, each with 
an opening angle significantly less than 180°. Allowing 



Formation of the double pulsar 



9 



TABLE 2 

Results from PSR J0737— 3039A Pulse Profile Geometry Fits 
FOR « < 90 ° . 



Pulse Height a 
(%) (°) 

(±68%) 


(5 Upper Limits 

n 

68% 95% 99% 


Pi 
(°) 
(±68%) 


Pi 

n 

(±68%) 


One-pole emission 


30 


90.31 


-11.5 
-11.3 


3.9 


26 


50 


115.8+0-* 




35 


90.21 


-17.1 
-18.2 


8.4 


55 


62 


115.0+O-t 




40 


90.41 


-16.8 
-16.7 


6.5 


38 


59 


114.1+0-? 




45 


90.41 


-18.2 
-18.2 


6.4 


33 


63 


113.0+0-? 




50 


90.31 


-16.2 
■16.3 


7.2 


40 


68 


111.91°;^ 




Two-pole emission 


30 


90.2111-3 


0.54 


3.9 


5.6 


-i-o-a_p.42 


12 l+'^-i 


35 


90.3l«8i 


2.3 


5.0 


5.9 


1 7 Q+3.5 
l'-«-0.37 


10.915J 


40 


90.21 


-16.3 
-16.2 


0.56 


3.5 


5.1 


17 s+5-8 
l'-°-0.97 


n 1+®" 

^J^-^-1.3 


45 


90.21 


-11.8 
-11.9 


0.43 


1.7 


3.1 


17 O+'l'^ 
^'■'^-0.87 


1 O"'"'^ '' 


50 


90.21 


-10.8 
-10.9 


0.47 


1.9 


3.6 


16 5+'*'^ 





emission to come from both magnetic poles in this pul- 
sar thus avoids the need for exotically large beams, or a 
reinterpretation of the the location of the emission beam 
center. 

It is for this reason that we strongly favor geometry 
such that this pulsar is a near- orthogonal rotator emitting 
from both magnetic poles, and with the spin axis closely 
aligned with the orbital angular momentum vector. From 
this point, we will thus use the results from the two-pole 
case in our discussion. We quote an average value for 
the spin- magnetic axis separation of a = 90±11±5°, 
where the first uncertainty quantifies statistical effects, 
calculated as the median value of the determined a un- 
certainties. The latter represents that which is due to 
systematic error (likely due to RFI and calibration effects 
as described in Section [5TT|) . and is taken to be the differ- 
ence between the maximum absolute uncertainty and the 
statistical error. We also quote 68%, 95%, and 99% aver- 
age upper limits to the pulsar A spin-orbit misalignment 
angle S of 0.85 °, 3.2 °, and 4.7°, respectively. 

The epoch of zero precession phase Ti is difficult to 
constrain, given the geometry of the system. It seems 
to favor values such that (j>so = ° or 180 ° falls within 
the span of our data set, denoted in Figure |4] by vertical 
dashed lines. However, for (5 ~ 0°, any fit value for Ti 
would be fairly meaningless, since it would be impossible 
to define a precession phase at any given time. Unless the 
geometrical parameters are very finely tuned, if pulsar 
A had a measurable misalignment angle S, we would ex- 
pect to see at least some noticeable evolution in the pulse 
profile shape in the 7 years (> 8% of its precession pe- 
riod) we have been monitoring this pulsar with the GET. 
This is an especially valid argument, given the significant 
changes one sees in other DNS profiles whose precession 
periods are far longer than that of PSR J0737— 3039A. It 
is for these reasons that we do not include Ti in Table [21 

Figure [5] summarizes these results for the z < 90 ° case. 
One can see that the results for all the fits we perform are 
consistent across all width data sets. Finally, regardless 
of which of the one- or two-pole emission models one 
prefers, we can see that the favored geometry for this 
system is one with a pulsar A spin axis at near right 



angle to its magnetic axis, with a very small spin-orbit 
misalignment. 

7. SPECIAL RELATIVISTIC ABERRATION 

iRafikov fc Lail ()2006l ) also suggested that differential 
delays due to latitudinal aberration should modify the 
pulse profile of pulsar A on an orbital time scale (see also 
iDamour fc Tavlodll992l : iStairs et aI]|2004D . This aber- 
ration arises from the fact that the line of sight for ob- 
servers in different inertial frames intercepts the emission 
beam at different co-latitude angles. Hence, the emission 
components located on either side of a hollow beam are 
expected to move towards and away from each other as 
the sight line intercepts the beam closer to the edge or the 
center, respectively. The orbital motion should therefore 
induce a periodic change of the latitudinal aberration, 
which will cause a sinusoidal variation in the separation 
between the pulse components. Thus, a given compo- 
nent of the emission beam will suffer a shift having an 
amplitude proportional to the orbital velocity multiplied 
by a geometric factor which depends on the beam and 
spin axis geometry: 



(3 



sin ( tan xo 



[cos i sin 77(cos ^ + e cos uj) 

— cosry(sin ^ -f- e sin w)], (9) 

where /3 is the orbital velocity, C, and rj have been defined 
in Section 15.21 ip is the orbital phase, e is the orbital 
eccentricity, lo is the longitude of periastron, and 



tan Xo 



sin a sin $ 







(10) 



cos a sin ^ — cos $o sin a cos C 

where in this case $0: defined as in Equation [1] is the 
spin longitude as measured from the beam center, and 
the other quantities here are defined earlier in this paper. 

If we believe that the two symmetric components of 
pulsar A's profile are the leading and trailing edges of a 
single pole, then they will suffer a periodic shift 2A<I>o- 
In the case of a two-pole model, one can show that under 
a symmetry assumption of the emission beams the shift 
will also be the same, though the definition of the beam 
geometry (e.g. a, p, etc.) will apply to the two poles. 
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From this, it follows that one can constrain pulsar A's 
geometry by measuring the separation between the two 
pulse components as a function of orbital phase. We have 
performed this analysis using obser vations at the GET , 
both with data from the SPIGOT (jKaplan et al.ll2005f ) 
and GASP backends. As is the case for the long-term 
evolution of the pulse profile, we found no significant 
change in the separation of the pulse components. We 
can set an upper limit to the value of the denominator 
of Equation [9l (sin^tanxo)"^, of 0.13, 1.33 and 4.68, 
corresponding to the la, 2a and 3(T confidence levels, 
respectively. Given that the equations for this analysis 
are similar to those found in Section 15.21 for the long- 
term profile evolution, the implications for the geometry 
of pulsar A are similar to those presented earlier, though 
less stringent. 

8. DISCUSSION 

These results are somewhat at odds with previ- 
ous work. Thr o ugh p osition angle variation analysis, 
iDemorest et al.l ()2004[ ) prefer solutions with a ^ 5° 
or a ~ 90°. They choose the former low-a solu- 
tion, in order to avoid an emission beam opening an- 
gle of approxima. t ely 18 0°. Although the analysis by 
[Manchester et all ()2005[ ) favors S < 60 °, consistent with 
our findings, they believe that a spin axis aligned with 
the orbital angular momentum is unlikely. This was due 
to early scintillation-based measurements of the trans- 
verse velocity of PS R J0737— 3039A/B, which supported 
a lar g e natal kick (iRansom et ahl l2004t iWillems et al.l 
120041: iColes et all l2005f ). As discussed earlier, how- 
ever, more recent and reliable timing measurements con- 
firm a transverse velocit y for the system o f less than 
lOkms"^. Furthermore. [Manchester et al.l (|2005[ ) arrive 
at best-fit values for a that virtually exclude a ~ 90°. 
This is in direct contrast to our findings, which strongly 
disfavor low-Q solutions (as well as those near 180 °) . 
As with lDemorest et al.l ()2004[ ). [Manchester et al.l(|2005[ ) 
enforce a half-beam opening angle of less than 90 ° in 
their model of the pulse width data, reasoning that val- 
ues larger than 90 ° are unphysical, or else means that 
the beam center has been misidentified. We have shown, 
however, that such a large beam can easily be avoided 
by allowing the emission to originate from both magnetic 
poles. 

Although we strongly prefer a two-pole model for po- 
lar cap emission from pulsar A, it has also been postu- 
lated that the observed wide pulse profile in many re- 
cycled pulsars is due to emission that originates in the 
outer magnetosphere. Here, caustic effects produce a fan- 
shaped beam that, depending on its orientation, could 
be vie wed with an ope ning angle that is greater than 
180° (|Ravi et al.l l2010t ). The observed coincidence of 
millisecond pulsar pulse profiles with counterparts found 
in observations with the Fermi Large Area Telescope 
(LAT) seeras to support this idea for ce rtain cases (e.g., 
lAbdo et al.ll2010l: lEspinoza et al.|[20T2l ). Such a model 
may also explain the frequent disagreement of recycled 
pulsa r polarization d ata with a standard RVM descrip- 
tion (|Yan et al.ll20lH ): for pulsar A in particular, it pro- 
vides reasoning for the observed symmetry in its pulse 
profile. However, analysis of Fermi LAT observations of 
PSR J0737-3039A/B clearly shows that the 7-ray and 
1.4-GHz radio data are not aligned in rotational phase. 



This, along with geometric modeling that uses both 7- 
ray and radio data, support a two-pole emission model 
(L. Guillemot et al. 2013, in preparation). Regardless of 
the choice of emission model, our main result of nearly 
aligned spin and orbital angular momentum axes for pul- 
sar A remains unchanged. However, the outer magneto- 
sphere emission model described above would give more 
relaxed constraints on the system geometry, similar to 
the one-pole emission model results shown in Table [2j 

A low spin-orbit misalignment angle means that a very 
smaU kick was likely given to the PSR J0737-3039A/B 
system due to the SN that resulted in pulsar B. This is 
supported by measurements of a low transverse velocity 
and small eccentricity for this system. The low mass of 
pul sar B al so lends crede nc e to the studi e s perf ormed 
bv IPiran fc Shaviv (12001 . IWillems et all pOOl . and 
iStairs et al.l (|2006i) . which favor a low pre-SN mass (< 
2 Mq) progenitor, and a low natal kick (< lOOkms"^) 
given the constraints of a low radial velocity. In the case 
of the latter study, our constraints are also consistent 
with their estimate of the misalignment angle, which they 
predict to be 0.5 < (5 < 11 ° (95.4% confidence). 

This analysis thus supports a scenario in which the 
PSR J0737— 3039B progenitor underwent a low mass- 
loss, relatively symmetric SN event. The prominent can- 
didates for such an event are a n ECS, or possibly the col- 
lapse of a low-mass iron core (iPods iadlowski et al.l[200l 
|2004), and we therefore favor one of these scenarios as the 
one that produced the double pulsar system as we now 
observe it. Discovery and observations of an increasing 
number of DNS binary systems will help to determine 
the prevalence of this type of DNS evolutionary channel. 

In addition, these new results provide an impetus for 
updating various studies that aim to constrain the kick 
velocity given to the PSR J0737-3039A/B system af- 
ter the second SN. For example, the simulations done by 
IStairs et al.l ([2006I ) can now be improved upon, by in- 
cluding our upper limits on i5 as a prior distribution i n 
their analysis. In the work done bv lWong et al.l (|2010( ). 
tighter constraints can also be achieved by using our S 
values, particularly if one adopts the two-pole emission 
model, as we believe is correct. The method of the latter 
study also holds the potential for constraining the radial 
velocity of this system. This is crucial if we are to once 
and for all determine whether the PSR J0737-3039A/B 
system velocity has a large radial component, the decid- 
ing factor in distinguishing between space velocities, and 
thus the magnitude of the natal kick on the system after 
the second SN. 

These results may also help to determine the cause of 
the very large misalignment angle for pulsar B reported 
by Breton et al. (2008), and more broadly, calls into 
question the origin of the pulsar B spin. For example, 
noting that the majority of the currently obse rved pulsar 
B spin must have come from the second SN, iFarr et aD 
(|2011h has addressed this issue by postulating that any 
natal kick resulting from the SN that formed pulsar B 
must have been displaced from the center of mass of the 
pulsar B progenitor star. The measurement of the spin- 
orbit misalignment in pulsar A presented here may con- 
tribute to tighter constraints on the magnitude of this 
displacement in such a model. 

Finally, we speculate on the possibility of other DNS 
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Fig. 5. — Summary of results of J0737— 3039A profile analysis 
assuming two-pole emission, for i < 90 ° . Each set of results is 
plotted against the corresponding fractional pulse height for which 
they were found. Top panel: mean pulse width. Uncertainties here 
are smaller than the plotting symbol sizes. Second panel from top: 
median value for a, with Icr uncertainties. Third panel from top: 
misalignment angle S. Here, the darkest to lightest regions repre- 
sent the 68%, 95%, and 99% upper limits, respectively. Bottom 
panel: median values of p, with 68% uncertainties. 



systems having been formed through an ECS (or sym- 
metric SN event) of the second-formed NS progenitors in 
these systems. If we look back at Table [Tl we see distinct 
similarities between several measured properties of PSRs 
J1756— 2251 and those of the double pulsar system. This 
is a close DNS systems, with a low eccentricity and a 
low second-formed NS mas s . It has been arg ued (e.g., 
Ivan den Heiiveil [200l [20071: lWong"eFaII [20To[ ) that the 
evolutions of this, and possibly other DNS systems, may 
also have undergone ECS events in their histories. 

PSR J1756— 2251 is the recycled, first-formed pulsar 
in that system. Long-term pulse profile monitoring ca n 
thus be done to constrain its geometry (jFerdmanl I2OO80 . 
Due to its proximity to the ecliptic, the determination 
of its transverse velocity has thus far presented a chal- 
lenge, however, the measurement of the right ascension 



component of its proper motion corresponds to a veloc- 
ity of 9.3 km s^^ in that directio n, based on a DM dis- 
tance of 2.5 kp c for this pulsar ( [Cordes fc Laziol 120021 : 
iFaulkner et"all l2005: Fcrdman 200l). This low measured 
velocity presents another possible clue in favor of the 
PSR J1756— 2251 system having undergone a similar evo- 
lution to the double pulsar. It is hoped that future ob- 
servations will elucidate whether this is indeed the case. 

It may be that ECS (or another small-kick SN scenario) 
may be the dominant channel for close DNS formation 
(jPodsiadlowski ct al. 2004) — the likelihood of a system 
remaining bound after two SNe is more likely if one of 
these involves a low mass-loss, small-kick event. On the 
other hand, the narrower range of the progenitor masses 
required for this channel to proceed may reduce the pop- 
ulation of suitable systems for this scenario. 

9. CONCLUSIONS 

In this article, we have presented the calculation of 
constraints on the geometry of the PSR J0737-3039A/B 
system through measurement of high signal-to-noise inte- 
grated profiles of pulsar A. Our results favor an interpre- 
tation where pulsar A's magnetic axis is nearly perpen- 
dicular to the rotation axis and where the system emits 
from both magnetic poles as a near-orthogonal rotator. 
Most interestingly, we find that the angle S between the 
pulsar A spin axis and the total angular momentum of 
the system is small, consistent with an alignment be- 
tween the two axes. This provides very important evi- 
dence in favor of pulsar B having formed as a result of a 
symmetric event, possibly an ECS. 

Further discovery and observation of other DNS sys- 
tems will help to determine how frequent such an evo- 
lution may be as a channel for DNS system forma- 
tion. Gaining a better understanding of the number 
of DNS systems that proceed through different forma- 
tion schemes will be important for population synthe- 
sis analysis, and for predictions of possible gravitational 
wave (GW) so urces for upc oming detectors such as Ad- 
vanced LIGO ( Harry et aLl ^2010) and Advanced Virgo 
(jAccadia et al.l 20 111) , which will be sensitive to GWs 
from DNS mergers. 

A final qualitative conclusion to be drawn from this 
work is that the remarkable long-term stability of the 
PSR J0737— 3039A pulse profile indicates that current 
and future timing analysis performed through observa- 
tions of this pulsar is likely to be extremely reliable over 
a long time span. This is especially important for new 
and continued tests of GR that will be undertaken with 
this astonishing system. 
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